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Abstract 

We consider initial-boundary problems for general linear first-order strictly hyperbolic 
systems with local or nonlocal nonlinear boundary conditions. While boundary data are 
supposed to be smooth, initial conditions can contain distributions of any order of singularity. 
It is known that such problems have a unique continuous solution if the initial data are 
continuous. In the case of strongly singular initial data we prove the existence of a (unique) 
delta wave solution. In both cases, we say that a solution is smoothing if it eventually becomes 
fc-times continuously differentiable for each k. Our main result is a criterion allowing us to 
determine whether or not the solution is smoothing. In particular, we prove a rather general 
smoothingness result in the case of classical boundary conditions. 

1 Introduction 

Solutions to hyperbolic PDEs demonstrate a wide spectrum of regularity behavior. The appear- 
ance of singularities in nonlinear cases is known as the blow up of a solution [1, 2]. Singularities 
can appear in a finite time even for small and smooth initial data [3]. In some cases, both 
linear and nonlinear, a solution with time either preserves the same regularity as it has on the 
boundary or becomes less or more regular in time. The singularities encountered in the latter 
case are called anomalous [4, 5]. Criteria for the appearance of anomalous singularities are given 
in [6, 7, 8, 5, 9]. These papers are devoted to the case of interaction between singularities weaker 
than the Dirac measure. The singularities resulting from this interaction turn out to be weaker 
than the incoming singularities. A different effect is observed in [11, 12], where the incoming sin- 
gularities are derivatives of the Dirac measure. In this case the interaction produces singularities 
stronger than the initial ones. We will focus on the phenomenon of improving regularity in the 
case of initial-boundary value problems with nonlinear local and nonlocal boundary conditions 
for first-order linear strictly hyperbolic systems. 

Specifically, in the domain IT = {(x, t) | < j; < 1, t > 0} we address the problem 

idt + A{x,t)d^ + A{x,t))u = g{x,t), {x,t) (1) 
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u{x, 0) = ifix), X G (0, 1) 

Ui{l,t) = hi{t,v{t)), l<i<k, te{0,oo) 
Ui{0,t) = hi{t,v{t)), k < i < n, tG(0, oo), 



(3) 



(2) 




(4) 



Note that boundary conditions (3) cover the cases of classical boundary conditions (if hi do not 
depend on v) and reflection boundary conditions of local and nonlocal type. 

We assume that 



for all (x, t) G 11. Condition (5) occurs in many applications, where the functions Uj for j < k 
(resp. k + 1 < j < n) describe the "species" that travel to the left (resp. to the right) and are 
reflected in x = (resp. x = 1) according to the boundary conditions (3). 

We will impose the following smoothness assumptions on the initial data: The entries of A, 
A, g, and h = {hi, . . . , hn) are smooth in all their arguments in the respective domains, while 
the entries of (p are allowed to be either continuous functions or strongly singular distributions. 

In the case of a continuous f (considered in Section 2), by a solution to problem (l)-(3) 
we mean a continuous solution, i.e., a continuous vector- function in 11 satisfying an integral 
system equivalent to (l)-(3). The existence and uniqueness of a continuous solution is proved 
in [13] (see Theorem 2). In the case of a strongly singular ip (considered in Section 4) by a 
solution we mean a delta wave solution, i.e., a weak limit of solutions to the original problem 
with regularized initial data that does not depend on a particular regularization. We refer the 
reader to [10, 9] for a more detailed definition and motivation of delta waves. Theorem 21 in 
Section 4 establishes the existence of a delta wave solution for a version of (l)-(3). 

It is clear that the regularity of initial conditions (2) constraints the regularity of a solution 
if the latter is considered in the entire domain 11. However, the influence of the initial data can 
be suppressed if the regularity behavior is considered in dynamics, starting from a point of time 
T. 

Definition 1 A solution u to problem (l)-(3) is called smoothing if, whatever m € N, there 
exists T > such that u e C"" {U n {t > T})"' . 

Our main result is a smoothingness criterion for solutions to problem (l)-(3) in terms of 
the layout of characteristic curves (Theorems 12 and 22). In the case of classical boundary 
condition, the criterion implies that the solution is smoothing whenever inf |Ai| > for all i <n. 
As another consequence, we obtain a class of boundary conditions under which the wave equation 
has smoothing solutions (see [14] for a special case of this result). 

Our analysis of problem (l)-(3) shows a phenomenon usually observed in the situations 
when solutions to hyperbolic PDEs change their regularity: the smoothness changes jump-like 
rather than gradually. Another feature of problem (l)-(3) shown in [13] is that, even if we 
allow non-Lipschitz nonlinearities in (3), this system demonstrates almost linear behavior. The 



Ai < • • • < Afe < < Afe+i < • • • < A, 



(5) 
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smoothingness effect in tfie non-Lipscfiitz setting contrasts witfi blowups [1, 2] observed in many 
nonlinear systems. 

In [6, 7], results similar to ours are obtained in some more special cases, namely for homoge- 
neous linear boundary conditions of local type with constant coefficients and continuous initial 
data. Some restrictions on (l)-(3) are imposed in [6, 7] by technical reasons, as the authors 
use an approach based on the Laplace transformation and the Green's function method. We 
extend these results to the case of general linear first-order hyperbolic systems, nonlinear nonlo- 
cal boundary conditions, and distributional initial data of the Dirac delta type and derivatives 
thereof. Note that we use a different approach based on the classical method of characteristics. 
This method suits well for understanding of the mechanism of the smoothingness effect. 

An essential technical difficulty to overcome in demonstrating the regularity self-improvement 
is caused by the fact that the domain of influence of initial conditions (2) is in general infinite. 
In other words, the regularity of solutions all the time depends on the regularity of the initial 
data. However, this dependence is different for the boundary and the integral parts of the 
equivalent integral form of problem (l)-(3). Since the boundary summands are compositions 
of boundary data with functions defining characteristic curves and hence are "responsible" for 
propagation of singularities, the smoothingness effect is encountered whenever the boundary 
summands have a "bounded memory" or, more rigorously, all characteristics of (1) are bounded 
and each boundary singularity expands inside 11 along a finite number of characteristic curves.^ 
Our strategy of obtaining the smoothingness criterion consists in identifying an appropriate 
class of boundary conditions ensuring the "bounded memory" property and in showing that the 
integral summands do eventually improve upon the regularity of the initial data. It is important 
for our analysis that the lower order terms, with the exception of the diagonal ones, contribute 
into the integral summands transversely to characteristic directions. This is ensured by (5). 
Therefore, the integral part of the system causes no propagation of singularities and, moreover, 
suppresses it. 

The mathematical motivation of the paper is the scope of the stability theory, the Hopf 
bifurcation analysis, and the investigation of small periodic forcing of stationary solutions of 
hyperbolic PDEs. The main reason why those techniques are well established for nonlinear 
ODEs and parabolic PDEs, but not for nonlinear hyperbolic PDEs, is that, in contrast to 
parabolic case, hyperbolic operators in general do not improve the regularity of their solutions 
in time (the question is closely related to propagation of singularities along characteristic curves) . 
This complicates, in particular, proving the Predholmness property of the linearizations which 
is crucial for the analysis of solutions to nonlinear problems. We provide a range of boundary 
conditions ensuring the desired smoothingness effect, which still makes possible to handle the 
bifurcation analysis of a class of nonlinear problems (this idea goes back to [15]). 

The practical motivation is caused by applications to mathematical biology [16], chemical 
kinetics (describing mass transition in terms of convcctivc diffusion and chemical reaction and 
analysis of chemical processes in countcrflow chemical reactors [17, 18, 19, 20]). and semicon- 
ductor laser dynamics (describing the appearance of self-pulsations of lasers and modulation of 
stationary laser states by time periodic electric pumping [21, 22, 23]). 

^This sharply contrasts with the case of a Cauchy problem where solutions cannot be smoothing because the 
boundary part all the time "remembers" the regularity of the initial data. 
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2 Continuous initial data 



Here we consider the case of continuous initial data <^{x). We will assume the zero-order com- 
patibility conditions between (2) and (3), namely 

<fi{0) = hi{0,v{0)), k + l<i<n, 

(^i(l) = hi{0,v{0)), l<i<k, ^ ' 

where v{Q) = ^(/?i(0), . . . , (fk{0),<fk+i{^), ■ ■ ■ ■, Vn{^)^ ■ By || • || we denote the Euclidian norm in 

Theorem 2 ([13, Thm 3.1]) Assume that the data Aj, Uij, gi, cpi, and hi are continuous func- 
tions in all their arguments, and the coefficients Aj are Lipschitz inx ^ [0, 1] locally int € [0, oo). 
Suppose that hi{t,z) are continuously differentiable in z & M" and for each T > there exists 
C > such that 

\\VMt,z)\\<C{\og\ogH{t,\\z\\))^/\ (7) 

where H is a polynomial in \\z\\ with coefficients in C[0,T]. If the zero-order compatibility 
conditions (6) are fulfilled, then problem, (l)-(3) has a unique continuous solution in 11 which 
can be found by the sequential approximation method. 

We now introduce the notions of an Expansion Path and an Influence Path, that will be our 

main technical tools. Let uJi{T;x,t) denote the characteristic of the i-th equation of (1) passing 
through {x,t) € IT. Let x be a characteristic of the i-th. equation of system (1). Suppose that x 
reaches dll at two points. Let (x, t) be that of these points having larger ordinate (hence, x = 
or a; = 1). We say that x reflects at {x,t) if 

dzhj{t, vi{t), . . . , Vi-i{t), z, Vi+i{t), . . . , Vn{t)) 7^ for some j < n and z G M. (8) 

In this case that of the characteristics Ljj{T;x,t) and ujj{T;l — x,t) which lies above the line 
{r = t} is called a reflection of x- If the reflection is defined by ujj{T; 1 — x,t), it will be called 
a jumping reflection. 

Remark 3 Note that condition (8) means that the i-th components of the vector v{t) participates 
in evaluation of Uj{0,t) for k -\- 1 < j < n and of Uj{l,t) for j < k. Whenever the continuous 
solution to problem (l)-(3) is known (see Theorem 2), condition (8) is easily checkable. Oth- 
erwise, one can use the following constructive sufficient condition: (8) holds true whenever for 
every (zi, . . . , Zj-i, Zj+i, . . . , Zn) £ W^~^ there is j < n such that dzihj{t, zi,. . . , Zn) ^ 0. 

Definition 4 A sequence of characteristics Xi)X2, ■ ■ ■ is called an Extension Path (EP) if each 
Xl+i is a reflection of xi- 

Definition 5 A sequence of curved segments Bi, . . . ,Bs is called an Influence Path (IP) if the 
following conditions are met. 

• Each Bi is a continuous part of a characteristic xi of the ji-th. equation for some ji < n; 
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Figure 1: An EP. 



• The whole path is monotone in the sense that the coordinate t continuously increases while 
moving along it; 

• The transition from Bi to .Bj+i can be of three types: 

- Bi and -Bj+i meet at a point {x,t) such that aji^^j^{x,t) ^ in any neighborhood of 

- Bi and -B^+i meet at a point {x, t) with x = or x = 1 and in any neighborhood of i a 
characteristic of the ji-ih. equation reflects to a characteristic of the j/+i-th equation; 

- Bi terminates at a point {x, t) with x = or x = 1, starts at the point (1 — x, t) 
on the opposite side of 911, and in any neighborhood of t a characteristic of the ji-th 
equation makes a jumping reflection to a characteristic of the j;+i-th equation. 

Roughly speaking, an IP is a piecewise continuous curve with smooth peaces Bi,. . . ,Bs lying on 
characteristic curves such that either Bi and meet within 11 or Bi terminates at {x,t) G 911 
and -Bi-i-i starts at the point (1 — x,t) on the opposite side of dlL. 

Remark 6 Note that any EP is an IP, while the opposite is not necessary true because segments 
of an IP not necessary lie on reflected characteristics. 

Definition 7 Define a set Xi, called the domain of influence of the initial data on Ui, as follows: 
{x,t) G Xi if the value Ui{x,t) can be changed by varying a function ip in (2). 

Since initial data expand inside 11 along characteristic curves according to boundary conditions 
(3) and the lower order terms of system (1), we have the following characterization. 

Lemma 8 (x, t) € Xi if and only if there is an IP emanating from the initial axis t = and 
going through a part of the characteristic LOi{T;x,t). 
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The sufficiency follows from the proof of the necessity in Theorem 12. The proof of the necessity 
is based on the constructive description of the domain of dependence of Ui at {x, t) which turns 
out to be the union of the IPs going through a part of uJi{T;x,t) below the line t = t. Under 
the domain of dependence of Ui at (xq, to) we mean the set of all points (x', t') G 11 such that if 
varying the function 

u{x, t') = tp{x) (9) 

in any sufficiently small neighborhood of x', the i-th component of the solution to problem (1), 
(9), (3) changes at {xo,to). 

Definition 9 Define a set X° to be the union of IPs emanating from the initial axis t = and 
satisfying Definition 5 with additional conditions imposed on the three transition types from Bi 
to Bi+i: 

- El and -B^+i meet at a point {x,t) such that Cj^^^j^ (x, t) ^ 0; 

- Bi and -B^+i meet at a point {x^t) with x = or x = 1 and at this point xi reflects to 
Xl+i'i 

- Bi terminates at a point {x,t) with x = or x = 1, starts at the point (1 — x,t) on 
the opposite side of 911, and at {x,t) the characteristic xi makes a jumping reflection to 
Xl+i- 

Note that Yf = X^. 

We now introduce two conditions that will occur in formulation of our results. 

(i) For every T > there exists T' > T such that, for all x € [0, 1], every EP passing through 
(x,T) lies below the line t = T' . 

(a) For every T > and i < n there exists T' > T such that, for all x S [0, 1] with 
(x,T) € X?, every EP containing the characteristic u)i{T;x,T) lies below the line t = T'. 

Remark 10 In many important cases conditions (i) and (a) can easily he reformulated and 
verified in terms of X, A, and h. Sometimes they can be verified even directly (see examples in 

Section 3). 

Set n'^ = {(x, t)\0<x<l,0<t<T}. 

Remark 11 Let T > 0. In the domain n\n^ let us consider problem (1), (9), (3) with aij = 

for all j (i.e., system (1) is decoupled) and with t' replaced by T in (9). Then condition (l) 
means that, whatever T > and ipix), the function ip{x) has a bounded domain of influence on 
Ui for every i < n. In other words, for any decoupled system (1), ifipix) is singular at some point 
X G [0, 1], then this singularity expands outside along a finite number of characteristic curves 
within n^' for some T' > T that does not depend on x (z [0,1]. In contrast to (l), condition 
(ll) means that, whatever T > and ^i{x), the function ipi{x) restricted to X? H {t = T} has a 
bounded domain of influence on Ui for every i < n. 

Theorem 12 Assume that the data Xi, aij, gi, and hi are smooth functions in all their argu- 
ments and conditions (5) and (7) are fulfilled. 
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• Sufficiency. Assume that condition (l) is fulfilled. Then the continuous solution to problem 
(l)-(3) is smoothing for any ip G C[0, !]"■ satisfying equalities (6). 

• Necessity. Assume that the continuous solution to problem (l)-(3) is smoothing for any 
if G C[0, 1]"^ satisfying equalities (6). Then condition (u) is fulfilled. 

Remark 13 One can easily see that the sufficient condition (i) implies the necessary condition 
(u), while the converse is not true. Nevertheless, in a quite general situation when X° = 11 for 
every i < n, it is not difficult to observe that conditions (i) and (ll) are equivalent. Hence we 
have the following result. 

Corollary 14 Assume that the data Aj, Oij, gi, and hi are smooth functions in all their argu- 
ments, X° = n for every i < n, and conditions (5) and (7) are fulfilled. Then the continuous 
solution to problem (l)-(3) is smoothing for any ip G C[0, 1]" if and only if condition (i) is 
fulfilled. 

Proof. Sufficiency. Assume that condition (i) is fuffiUed. Define a sequence Ti,T2,... 
inductively by the following rule. Let Ti be the infimum of those r > for which there is x and 
an EP passing through (x, 0) and lying below the line t = t; let Tj for j > 1 be the infimum 
of those r > for which there is x and an EP passing through (x,Tj_i) and lying below the 
line t = Tj-i. Note that Tj is monotone and approaches the infinity. The latter fact is a simple 
consequence of the smoothness assumptions on A. By Theorem 2, problem (l)-(3) has a unique 
continuous solution u in n. It suffices to show that u G (ll \ 11^^ )". 

Consider problem (l)-(3) first in 11 \ II^i. The solution satisfies the system of integral 
equations 



Ui{x,t) = hi{ti{x,t),v{ti{x,t))) 



(10) 



+ 



j Ei(T;x,t) 



ti(x,t) 



dr, i < n, 



^=uJi{T;x,t) 



where Ei(T;x,t) = exp J^^ aji(cjj(Ti; .t, t), ri) dri and ti{x,t) denotes the smallest value of r > 
at which this characteristic reaches OH. In the sequel, along with the equation ^ = a;i(r; x, t) we 
will also use its inverse form t = oji{^;x,t). Due to the definition of Ti, in the boundary term 
hi{ti{x,t),v{ti{x,t))) in (10) we can substitute 



Vj{ti{x,t)) = hj{tj{xj,ti{x,t)),v{tj{xj,ti{x,t)))) 

ti{x,t) 



(11) 



+ J Ej{T;xj,ti{x,t)) 



■ XI ^i'^^' t)us{^, r) + gj{i, t) 



dr, 



i=ojj(T;xj,ti(x,t)) 



where Xj = Oifl<j<k and xj = 1 if k + 1 < j < n. Continuing in this fashion, the right-hand 
side of (10) can eventually be brought into a form depending neither on v nor on ip. This version 
of (10) will be referred to as (10'). Wc begin with establishing the (7^ (lI \ 11"^^ ) -smoothness of 
u. It will be proved once we show that the right-hand side of (10') has a continuous partial 
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X 



tJl,tj(Xi,ti(x,t))) 



(a.(v,0,tf(x,t)) 

J . 










ti(0,ti(x,t)) 




tj(Xi,t.(x.t)) 


/ (Ojx;l,tj(0,ti(x,t))) 





t 
Figure 2: The domain of integration. 



derivative in x. The latter can be done by transforming all integrals occurring in (10'). The 
transformation of each integral follows the same scheme, which we illustrate by example of the 
integral expression 



ti{x,t) 



Iijm{x,t)= J Ej{T;Xj,ti{x,t))ajm{^,T)Um{^,T) 



dr. 



(12) 



We will use assumption (5). Suppose, for instance, that i < k, j > k + 1, m < k, and 
tmi^,tj{xi,ti{x,t))) > 0. This entails, in particular, that Xi = and Xj = 1. Due to (10), we 
obtain 



ti{x,t) 



Iijm{x,t) = J Ej{T;0,ti{x,t))ajm{^,T) 

tj{OM{x,t)) 

T 

+ j £;^(ri;^,T)|-^( 



hm{tmi^,T),v{tm{^,T))) 



ti{x,t) 

^ / 



dri 



Z=U)m{Ti;C,T) 



dT 



j Ej{p; 0, ti{x, t))ajmiu;j{p; 0, ti{x, t)),p) 



p=9{0,T;x,t) 



hm{T,v{T))dT 
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+ / Ej{p;0,ti{x,t))ajmir], p)E„ 

" Sm. 



r]=u}m{6{C,T;x,t);^,T),p=e{^,T;x,t) 



X < 



(13) 



where Sm is the area shown in Fig. 2 and 9{£,,T;x,t) denotes the t-coordinate of the point 
where the characteristics ojj{Ti;0,ti{x,t)) and uJm{Ti;^,T) intersect. The other cases are sim- 
ilar. For example, if m > k + 1, then in the formula (13) index m should be replaced by I 
and the integration over [tm{^,tj{0,ti{x,t))),ti{x,t)] should be replaced by the integration over 
[tjiO,tiix,t)),tiiO,ti{x,t))]. 

The desired C^-smoothness of Iijm{x,t) follows from the C^-smoothness of 9{^,T;x,t). 
The latter is a consequence of the smoothness properties of A. Indeed, from the equality 
u}j(9{^,T;x,t);0,ti{x,t)) = LOm{G{^,T", x,t);(,,T), we conclude that, \i dx9 exists, then it is given 
by the formula 



dx9{^,T;x,t) = 



dojj {9{^, T;x,t);0, U {x,t)) dU {x, t) 



dti dx 
du)m{9{^, r; x, t);^, r) dojj{9{^, r; a;, t); 0, ti{x, t)) \ 



= —Xj{0,ti{x,t))exp < 



89 89 

' e{^,T:x,t) 

, ti{x,t) 
' ti{x,t) 



da 



uij{cr;0,ti {x,t)) 



X {Xi{Q,ti{x,t))) ^exp < 



\^(C,^) 



da 



€,=0Ji{a;x,i) 



>^m{^m{Oii, t; X, t);^, t),9{^, t; x, t)) 



1 -1 



- \j {ujj{9{^, r; X, t); 0, U {x, t)),9{^, r; x, t)) 

Thanks to the equality ujm{9{^,T;x,t);(^,T) = ujj{9{^,T;x,t);0,ti{x,t)) and condition (5), the 
function Iijm{x,t) is continuously diffcrcntiable in x. 

Thus, the right-hand side of (10') is continuously differentiable in x. Therefore, dxU G 
C (n \ n-^i)". The membership of u in now directly follows from system (1). 

In the next step, we prove that %i G C'^ tQ \ 11-^^)". For dxU wc have equations 



dxUi{x,t) = {Riu){x,t) + J Ei{T;x,t) 

ti{x,t) 

-d^XiiS,, T)d^Ui{^, t) + (d^fi) (^, r, u 



(14) 



dr. 



i < n, 



^=0Ji{T;x,t) 



where fi{^, r, u) = gi{^, r) - Yyj=i aij{^, r)uj, 

{Riu){x,t) = X['^{yi,T) fi{yi,T,u) -Vyhi{T,v) • w'(t) - {drhi){T,v) 



T=ti{x,t) 



(15) 
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s+3 

dr 



= for A; + 1 < z < n, and j/j = 1 for 1 < ? < A;. Here w ■ z denotes the scalar product in 
In (15) we can represent 'y'(i) in the form 

^K*) = /i(a^i>*>^) - Xj{xj,t)dxUj{xj,t) = fj{xj,t,u) - Xj{xj,t) {Rju){xj,t) 
t 

+ J Ej{T;xj,t)^-^ajs{^,T)d^Us{S,,T) - d^\j{^,T)d^Uj{^,T) 
+ id^fj){^,T,u)) 

We continue in this fashion up to getting a representation of the boundary term (15) that does 
not depend on v', what is possible due to the definition of T2. To show that the right-hand side of 

the obtained expression for dxUi{x,t), say (14'), is continuously differentiable in x, we transform 
all integrals contributing into (14') similarly to (13). Using the fact that u G C^(n\n-^i)", one 
can easily show that u G C^'^ (ll \ H^^)" The desired (ll \ n^2)_g]^QQ^];ijiggg Qf ^j^g solution 
is then a direct consequence of system (1) and its differentiations. 

We further proceed by induction on k. Assuming that u G C^~^ (ll \ n-^'^-i) for some k >2, 
we will prove that n G (7'= (n\n^'=)". We differentiate (1) k times in X, thereby obtaining a 
system for d^u. The integral form of this system is similar to (14). Analogously to (14), the 
definition of makes possible an integral representation of the system for d^u which does not 
depend on v^''^ and includes integrals of d^^^u similar to Iijs{x,t). To show the C^-smoothness 
of the solution, wc transform the integral terms analogously to (13). Finally, the C'^-smootlmcss 
of u outside of 11-^*= follows from suitable differentiations of system (1). The sufficiency is thereby 
proved. 

Necessity. Suppose that condition (ll) is not fulfilled and prove that the solution to problem 
(l)-(3) is not smoothing for some (p G C[0, 1]". Fix T = to > and i < n such that for all 
T' > to there is an EP containing the characteristic LOi{T]x,to) for some {x,to) G Xi and going 
beyond 11^'. 

Since all singularities of solutions expand along EPs, it is sufficient to prove that there exist 
m G N and ip G C[0, 1]" such that, whatever (xo,io) € X?, the solution u is not C""-smooth at 
{xo,to)- 

By Lemma 8, for any (x, to) G there is an Influence Path Bi, . . . ,Bs emanating from the 
line t = and going through a part of uJi{t;x,to)- Denote the smallest possible "length" of such 
a path by s{x). By the smoothness assumptions on the initial data, if x' is sufficiently close to 
X and (a;',to) £ X^., then s{x') < s{x). Since n {t = to} is closed, the standard compactness 
argument implies that s{x) is bounded by a constant sq uniformly over all x. 

Let be a continuous nowhere differentiable function. Consider an arbitrary {xo,to) G X°. 
Let us fix a shortest Influence Path Bi, . . . , Bg from some {zq,0) to (xo,to) with Bg lying on 
u!i{T; xo,to) and with the transition from Bi to -B/+1 as in Definition 9. Using the smoothness 
assumptions on A and h, we can suppose that < < 1. We now intend to prove that 
the solution u is not C®-smooth at (xo,to)- Since s < so, this will give us the necessary part 
of the theorem. Let (zi,Ti) denote the starting point of B^+i. Suppose that Bi is a part of 
a characteristic xi of the j'^-th equation. As ipj-^ is not continuously differentiable at zo, the 
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Figure 3: The set B. 

function Uj^ is not C^-smooth along xi due to the definition of a characteristic. If X2 is a 
reflection of xi, then Uj^ is not differentiable at (zi,Ti) and hence Uj^ is not differentiable along 
X2- Otherwise, by Definition 5, we have aj^j^{zi,Ti) / and Uj^ is not C^-smooth along X2-, 
because the integral of cij2ji%i along X2 in the integral form of the j'2-th differential equation is 
not a C^-function and this nonsmoothness cannot be compensated by any other summands in 
the integral representation of our problem. It follows that Uj^ is not C^-smooth, in particular, 
at {z2,T2). Continuing in this way, we arrive at the conclusion that Ui is not C^'-smooth along 
0Ji{t\XQ,tQ) and hence u is not C*-smooth at (xo,to)- 

Going into the details of the above argument, let us represent Ui{xQ,tQ) in an integral form 
with an integration over a neighborhood of the Influence Path Bi, . . . , Bg. We focus on the 
case of s = 3 where characteristic X3 is not a reflection of X2 and X2 is not a reflection of 
Xi (see Fig. 3). A similar or eve n simpler argument works as well for other possible cases. 
Extend each Bi to B'^ = Xi ^ Il^K Given e > 0, let B denote the union of B'g and one-sided 
neighborhoods of B[, . . . ,B'g_i, bounded by the characteristics Ljjj,^{T]UJi{T2 — e; XQ,to),T2 — e) 
and Ldj-^ (r; loj^ (Ti + e; uJi{T2 — e; xq, to); ?2 — e), + e). Fix an e > so that B contains neither 
(0,0) nor (1,0). 

Now we write an integral representation of Uj(xo,to) in terms of u over B. We start from 
the formula (10) for Ui{xQ,tQ) and rewrite a part of the integral in the right-hand side as 
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T2 



T2-e 



T2 



Ei{T; xo,to)aij2{C{r),T) 



T2-e 



hn {tjAiij),T),v{tjMr),r))) 



T 

+ j Ej^in;^{T),T) ^ aj2,(a;j2(ri;^(r),r),ri)u,(a;j2(ri;C(r),r),ri) 



dri 



dr. 



(16) 



By construction, aij2{z2,T2) ^ 0. Furthermore, we consider the part of the integral in (16) over 
the area denoted in Fig. 3 by S and transform it similarly to above as follows: 



T2 



Tl+6 



T2-e Ti 

X aj2ji {ujj^ (ti ; ^(t) , r) , Ti)uj^ (uj^ (n ; ^(r) , r) , n ) dn dr 

T2 Ti+e 



J Ei{T;xo,to)aij^{(,{T),T) J Sj2(ri;^(r),r)aj2j,(wj2(ri;^(r),r),ri) 

2-e Ti 

n 

jiKi(0;'^j2(n;C(T),r),ri)) + J Ej,{T2;ujj2{Ti;^{T),T),Ti) 



X] (^hsi^jl {r2;0Jj2iTi;^{T),T),Tl) ,T2)Us{u)j, {t2; OJj^in; ^{t) , t) ,Tl) ,T2) 



+9h{^h ('r2;wj2(n;C(T),r),ri) ,r2) 



dT2 



dri. 



Again by construction, aj^j^{zi,Ti) ^ 0. Combining (10), (16), and (17), we see that 

Ui{xO,to) = Qe[uKxo,to) + $f(xo,to), 

where is a certain operator and 

T2 Ti+e 

M{xo,to)= J Ei{T;xo,to)aij2{^iT),T) J Ej^{Ti;£,{T),T) 

T2-6 Ti 

Xfljwi {^j2iTi;^{T), r), ri)<y9ji (wji (0; u^j^ (n; ^(r), r), n)) dn dr. 



(17) 



(18) 
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A similar representation for Ui{x,t) holds in a sufficiently small neighborhood of {xo,to). From 
the derivation of (18) it follows that Ui at (a:o,to) cannot be more regular than any of the 
summands in (18). Hence it suffices to show that $|(a;o,to) is not smooth at {xo,to). Let us 
make simple transformations: 

T2 

^t{xO,tQ)= J Ei{T-XQ,tQ)aij^{^{T),T) 

T2-e 

viTi,r) 



X 



J EjiiPix, t); C(t), T)aj2ji {a{x, t),P{x, T))(pj^ (x) dx dr 



v{Ti+e,T) 

= J Ei{T;xo,to)aij^{^{T),T) 

T2-e 



J Ej2 W{x, T)-^{T),T)aj2j^ {a{x, t),P{x, T))(pj^ (x) dx 




+ J Ej^{l3{x,Ty,^{T),T)aj2j^{a{x,T),/3{x,T))(pj^{x)dx 



dr 



r]{Ti+£,T) 
n{Ti,T2) 

Ei{p{y);xo,to)aij^ {j{y),p{y)) 

ri{Ti,T2-e) 

y 

X 



»7(Ti+£,T2) 



j EjMx, P{y));7{y),p{y))aj2ji ((^i^, p{y)),Pix, p{y)))ipji (x) dx dy 



v{Ti+e,T2) 

+ j Ei{C{y);xo,to)aij^{a{y),C{y)) 



v(Ti+s,T2-e) 


X J EjMx, Ciy));<^iy), C,{y))o.j2h i^i^^ C{y)),Kx, C,{y)))vh i^) dxdy, 

y 

where rj{t,T) = ujjj^{0;ujj2{t;£,{T),T),t); (a(a;, r), /3(.t, r)) denotes the intersection point of the 
characteristics (ri; x, 0) and a;j2(Ti; ^(r), r); {'y{y), p{y)) denotes the intersection point of 
the characteristics Ui{T;xo,to) and u}j2{T]U)j-^^{Ti]y,0),Ti)] {cr{y),C{y)) denotes the intersec- 
tion point of the z-th and the j2-th characteristics passing through the points a;j(r;xo,to) and 
Uj^ (r; (Ti + e; y , 0) , Ti + e) . Note that a, /3, 7, p, a, and ^ are (at least) continuous functions, 
what easily follows from the smoothness assumptions imposed on the initial data. It follows 
that Ui is not C'^-smooth at {xo,to) (even if a, /3, 7, p, a, and ^ are C°°-functions). 

The necessity is proved. □ 
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3 Examples 



Here we give some examples to show how the criterion given by Theorem 12 works. Each of these 
examples is rather general and interesting by its own. Throughout this section it is supposed 
that all characteristics of system (1) are bounded. This assumption is not restrictive from the 
practical point of view. It is true, for example, whenever inf |Ai| > for all i < n. 

3.1 Classical boundary conditions 

As a partial case of problem (l)-(3), consider (1), (2) with classical boundary conditions 



Theorem 15 Assume that the data Xi, aij, gi, and hi are smooth in all their arguments. Sup- 
pose that condition (5) is fulfilled. Then the continuous solution to problem (1), (2), (19) is 
smoothing for any ip G C[0, 1]" satisfying (6). 

This result is a straightforward consequence of Theorem 12. Indeed, since no reflection from 
the boundary is possible, every EP passing through {x,t) consists of a single characteristic. 

Moreover, due to the smoothness assumptions on A and the definition of a characteristic, for 
any T > there is T' > T such that all characteristics passing through the line t = T \e\ below 
t = T'. Hence the sufficient condition (l) is fulfilled and the theorem follows. 

Remark 16 Even in the case of smooth classical boundary conditions (19), the domain of in- 
fluence of the initial data on Ui for every i < n in general is unbounded ( due to the lower order 
terms in (1)). In spite of this, the influence of the initial data on the regularity of u becomes 
weaker and weaker in time causing the smoothingness effect. 

3.2 Periodic boundary conditions 

Suppose that at least one component of the solution, say the first one, satisfies periodic boundary 
condition. Specifically, the first equation in (3) is written in the form ui (0, t) = ui{l,t). Then the 
domain of influence of initial data (2) on ui is the whole U. Moreover, for each {x,t) G H, there 
is an unbounded EP passing through (x, t) (at least one constructed by means of characteristics 
of the first equation). This entails that the necessary condition (ii) is not fulfilled and hence the 
solution to problem (l)-(3) is not smoothing. 

3.3 Nonsepeirable lineEir boundairy conditions of local type 

Consider now the reflection boundary conditions 



UiiO,t) = hi{t), k + l<i<n, t G (0, oo) 
Ui{l, t) = hi{t), 1 <i <k, t e (0, oo). 



(19) 



k 




k-\-l <i <n, i G (0, oo) 



(20) 



n 




t G (0,oo), 



j=k+l 
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where bij and Cij are constants. Systems like (1), (2), (20) cover linearizations of many mathe- 
matical models for chemical kinetics [19, 20]. 

Our aim is to reformulate sufficient condition (i) constructively in terms of data (specifically, 
in terms of boundary data) of our problem. Note first that, if bij ^ for some k + 1 < i < n 
and j < k (resp., ^ for some i < k and A; + 1 < j < n), then each characteristic of the 
j-th equation reflects from the boundary a; = (resp., from the boundary x = 1), the reflections 
being characteristics of the i-th equation. Since all EPs passing through {x,t) are constructed 
by means of subsequent reflections, for each EP passing through (x, t) and consisting of more 
than one smooth piece, there is a unique sequence (finite if the EP passing through {x, t) is 
bounded and infinite otherwise) either of kind 



or of kind 



Ci2«l ' ^*3j2! ''*4*3' ^isM! • • • (2^) 



with nonzero elements. Condition (t) is satisfied iff for every T > all EPs passing through 
{x,T) are bounded uniformly in x G [0,1]. The latter is equivalent to the statement that all 
the sequences (21) and (22) are finite. This means that condition (<,) is expressible here in the 
algebraic form, namely as a finite number of equalities of kind 

for all possible such that the entries of matrices Bn-k,k = {bij}k+i<i<n,j<n and 

Ck,n-k = {cij}i<k,k+i<j<n appear in (23) no more than once. 

Condition (23) can be easily reformulated in the following form (see also [7]). Set 



R 



Ck,n—k Ik 



where Ij is a unit matrix of size j. Consider an expansion of the determinant of R along the 
first n — k rows, namely, 

deti?= Ai,...,^„_,i^n,...,^„_,. (24) 

n<---<i„_fc 

Here the respective determinants of order n — k are denoted by Di^,. and their cofactors by 
Fii,...,i„-k- turns out that each summand in (24) starting from the second one is a product of 
kind (23). One can easily show that (23) is fulfilled iff all summands in (24) excepting for the 
first one, which is obviously equals to 1, vanish. Note that the latter condition is obtained in [7] 
for problem (1), (2), (20) with coefficients in (1) depending only on x. 

3.4 Smoothingness phenomenon for IBVPs for the wave equation 

Using Theorems 2 and 12, one can immediately obtain correctly posed IBVPs for the (nonho- 
mogeneous) wave equation 

{d^-a'dl)u = f{x,t) (25) 
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and state a smoothingness result for them. For instance, let us consider (25) subjected to initial 
conditions 

u{x,0) = (p{x), dtu{x,0) = ijj{x) (26) 

and boundary conditions 

u{0,t) = hi{t,u{l,t),{dtu + ad^u)\^=o) , 

{dtu + adxu)\x=i = h2{t,u{l,t),{dtu + ada:u)\x=o) . 

The problem (25)-(27) is equivalent to the following problem for the (2 x 2)-first-order hyperbolic 
system 

{dt + adx)u = w, {dt - adx)w = f{x, t) (28) 

u{x,0) = (p{x), u;(x, 0) = ^(x) + a^^^ (29) 
u{0,t) = hi{t,u{l,t),wiO,t)) , 

w{l,t) = h2{t,u{l,t),w{0,t)) . ^^^^ 

Under a continuous solution to problem (25)-(27) we mean the continuous solution to problem 
(28)-(30). In the next theorem conditions (l) and {u) are supposed to be formulated correspond- 
ingly to problem (28)-(30). The following smoothingness result for (25)-(27) is a straightforward 
consequence of Theorems 2 and 12. 

Theorem 17 Assume that the data f and hi are smooth functions in all their arguments and 

conditions (5) and (7) are fulfilled. 

• Sufficiency. // condition (i) is true, then the continuous solution to problem (25)-(27) 
is smoothing for any (p € C^[0, 1] and tp G C[0, 1] satisfying the zero-order compatibility 
conditions between (29) and (30). 

• Necessity. Assume that the continuous solution to problem (25)-(27) is smoothing for any 
ip G C^[0, 1] and ip G C[0, 1] satisfying the zero-order compatibility conditions between (29) 
and (30). Then condition (a) is fulfilled. 

4 Distributional initial data 

Now we address the case of distributional initial data. Specifically, we investigate system (1) 
with initial conditions 



Ui{x, 0) = bi{x) = bir{x) + bis{x) = bir{x) + ^ Cij5^^ii\x - x*^), 

i < n, < a;*i < • • • < a;*„. < 1, a; G (0, 1) 
and boundary conditions 

n 

Ui{(},t) = J2Pi3it>jit) + ri{t,v'^{t)), k-\-l<i<n, iG(0,oo) 

n 

Ui{l,t) = J2Pijit>jit) + ri{t,v''{t)), l<i<k, tG(0,oo), 



(31) 



(32) 
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where the regular part br{x) = {bir, ■ ■ ■ , hnr) of the initial data h{x) = (61, ... , 6„) is a contin- 
uous vector-function on [0, 1], ri{t,y) € ([0, 00) x R"), n and VyVi are bounded in y G 
uniformly in t over any compact subset of ]R_|_, and Pij{t) € C°°[0, 00). Also, Cjj G M, mj G N, 
and lij € Nq. Furthermore, v^{t) just redenotes the vector-function v{t) introduced in Section 1 
by (4). We will use such notation below with other superscripts. Without restriction of gen- 
erality we can suppose that the zero-order compatibility conditions between (31) and (32) are 
fulfilled, namely, 

n 

brM = J2Pij(^H(^) + r^{0,v''iO)), k + l<i<n 

(33) 

hri{l) = ^pi,(0)t;j(0) +ri(0,i;"(0)), l<i<K 

where v"(0) = (6r-i(0), . . . , &r-fe(0), 6r,fe+i(l), • • • , hmO-)). 

Note that system (1), (31), (32) is unsolvable in the classical sense, since nonlinearities in the 
boundary conditions and singularities in the initial data entail nonlinear operations in the space 
of distributions. We, therefore, aim at constructing a delta wave solution to (1), (31), (32). For 
this purpose we regularize the singular part hs{x) of initial data (31) by means of the convolution 
with the so-called model delta nets {ip^)^^Q where ip G 'C(M), J ipdx = 1, ipe{x) = l/e(p{x/£), 
and J \(ps\dx is bounded by a constant independent of e. The regularized problem for is this: 

{dt + Ad^ + A)u' = g{x,t) 

mi 

ul{x, 0) = bir{x) + h%{x) = hir{x) + ^ Cijip-t'\x -x*j), l<i<n 

n 

nf(0,t) = ^^^.,(t)^;;''(t) + r,(i,i;"'^(t)), k + l<i<n (^^) 

n 

<(l,t) = Y.MiK'it) + ri{t,v'''%t)), l<i<k. 
3=1 

Here v"'^(t) = (uf (0, t), . . . , u|(0, t), u|_^i(l, t), . . . , <(1, t)). Decompose A = D + F into the 
diagonal and off-diagonal parts D and F, respectively. Our aim is now to show that -^e^o 
z+w m 2?'(n) where a singular part z and a regular part w are solutions to systems, respectively, 

{dt + Ad^ + D)z = 

z{x,0) = bs{x) 

n 

ZiiO,t) = ^pij{t)vj{t), k + l<i<n 

i=i 

n 

Zi{l,t) = J]p,,(t)i;|(t), l<i<k 



17 



and 



{dt + Ad^ + A)w + Fz 
w{x, 0) 



g{x,t) 

br{x) 



n 



Wi{0,t) 



Y,Pij{t)vJ{t)+ri{t,v''), k + l<i<n 



(36) 



i=i 



n 



Wi{l,t) 



Y,Pij{t>f{t)+ri{t,v^), l<i<k 



(the two systems are called nonlinear splitting). System (35) is linear, corresponds to the singular 
part of the initial problem, and is responsible for propagation of singularities. System (36) is 
nonlinear and corresponds to the regular part of our problem. Note that, since F is the off- 
diagonal part of A, the singular forcing term Fz contributes into the hyperbolic system (36) 
transversely to the characteristic directions and, therefore, cannot produce strong singularities 
for w. We now show that problems (35) and (36) are uniquely solvable. 

Definition 18 By J^, we denote the union, over i < n and j < nii, of all EPs containing the 
characteristic coi{T;x*j,0). Furthermore, J will denote the union, over i < n and j < rui, of all 
EPs emanating from the point {x*j,0). 

Note that J* C J. 

Lemma 19 System (35) has a unique distributional solution z representable as the sum of 
singular distributions concentrated on characteristic curves in the set J* . 

The lemma can be easily proved by the method of characteristics. 

Let C (n \ j) be the space of piecewise continuous functions on 11 with (possibly) first-order 
discontinuities on J. 

Lemma 20 System (36) has a unique C (ll \ -solution w. 

Proof. The uniqueness can be easily proved by considering the problem for the difference 
of two C (n \ j) "-solutions to (36) and applying the method of characteristics. 

To prove existence, consider w = w + w, where iB is the solution to the linear system with 
the strongly singular forcing term 



{dt + Kd^ + D)w + Fz = g{x,t) 

U)(x, 0) = br{x) 



n 



(0,t) = ^pij{t)vf{t), k + l<i<n 



(37) 



n 



Wi{ 



(l,i) = J2pij{t)vf{t), l<i<k 
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and w is the solution to the nonlinear system 

{dt + Ad^ + A)w + Fw 
w{x, 0) 

Wi{l,t) 

Note that, if it; is a solution to (37) and w is a solution to (38), then w + w is a solution to (36) 
indeed. 

To show that system (37) has a (unique) solution w G C \ j) , we use the method of 
characteristics again, thereby representing w explicitly in the integral form. Since the singulari- 
ties in the forcing term Fz are strong and transverse to the corresponding characteristics of the 
hyperbolic system, the Volterra integrals of Fz in the integral representation of W are at worst 
discontinuous functions with the first order discontinuities (if any) on J. 

To show that problem (38) has a (unique) solution w E C \ j)" , we follow the proof of 
Theorem [13, Thm 2.1] with minor changes. Fix an arbitrary T > and rewrite (38) in XI-^ in an 
equivalent integral form with the boundary term not depending on v'^ (applying a finite number 
of integrations along characteristic curves up to reaching the initial axis) . Using the assumption 
that Vyri{t,y) is bounded locally in i G M_|_ and globally in y G R", it is not difficult to check 
that the operator defined by the right-hand side of the integral form of (38) maps C (ll \ j) " 
into itself and is contractible on II^o for some £ T. The local existence-uniqueness result in 
n^o then follows from the Banach fixed point theorem. If 6q < T, we consider problem (38) 
on \ n^o with the initial condition replaced by u{x,9o) = ip{x) G C ([0, 1] \ J)", where the 
function il;{x) is determined in the first step. Similarly to the above, the right-hand side of the 
integral form of this problem defines the operator mapping C ((ll \ 11^°) \ j) into itself. Set 
n'' = 0. Since the contraction property of the operator under consideration docs not depend on 
the initial conditions, one can choose the value of from the very beginning so small that the 
operator is contractible on n(*+i)^o \ 11*^° for any < s < \T/9o] — 1. Applying the Banach 
fixed point theorem with respect to the domain II^^o \ n^°, we conclude that problem (38) has 
a unique C (^(^tf ^" \ n'^o) \ ''-solution. The proof of the unique solvability in n^^+i)^" \ 
for each s G {2, . . . , IT/Oq] — 1} goes over the same argument. Since T > is arbitrary, the 
unique solvability of (38) in C (n \ j)" follows. 

We have thus proved that (36) has a unique solution w in C(n\ j)", which equals to 
H) + w. □ 

We are now prepared to state our result about delta waves. 

Theorem 21 Assume that the data Xi, aij, gi, pij, Vi, and bir are continuous functions in all 
their arguments, Xi are Lipschitz in x € [0, 1] locally in t £ [0, oo), and both ri and Vyri are in 
L°° ((0,r) X R") for any T > 0. Given e > 0, let , z^ , and be the continuous solutions 
to problems, respectively (34), (35), and (36) with z^ , , and 6f in place of z, w, and hg, 
respectively. Then 








n 



n 



(38) 



Y.Pij{t)vf{t) + ri{t,v'' + v^), \<i<k. 
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1. - -w^ in Lj^cO^) as e ^ 0; 

2. ^ z in S'{U) as £^ 0; 

3. ^ w in C{K) os e — > for any compact subset K C Il \ J; 

4. ^ z + w in X>'(n) as £ — >■ 0, 

where z G £"'(11)" and w E C are the solutions to problems (35) and (36), respectively. 

Proof. 1. Let Jl be the union, over i < n, j < nii, and a G [— e;e], of all EPs containing 
oji{T]x*j + a,0). Note that the set consists of the "tubes" forming the support of z^ . Set 
= Jl n {dli \{t = 0}) . 

To show that the splitting of into a "regular" and a "singTilar" parts and z^ is correct 
(what is stated in Item 1 of the theorem), we consider the following problem for = —z^ —w^: 

{dt + Kd^ + A)a^ = 
a=(x,0) = 

n 

af(0,t) = ^pi,(t)i;f'^(i) + r-,(i,i;"'^)-ri(i,i;"'''^), k + l<i<n (39) 

n 



Rewrite 



ri(t,7;"'^)-r,(t,i;-'^) 



1 



j {Vyn) {t, av''+^+''' + (1 - a)v^+''') da ■ v""'' + {nit, v'"+''') - n{t, v^'')) 





n 



^Y.^mvf'+m), (40) 

where the symbol = introduces the short-hand notation Kfj and L|. By assumption, rj(t, y) and 

yyri{t,y) are in L°° ((0,T) x M") for any T > 0. Hence Kf-{t) and -Z^f(i) are bounded uniformly 

in t G [0, r] and e > for an arbitrarily fixed T > 0. 

To prove the desired statement, wc will use dominated convergence. It is sufficient to show 

2^ 

that, given T > 0, the sequence is bounded on 11 uniformly in s: > and a^{x,t) — )• as 
£ — >^ pointwise off J*. 

J' 

Fix an arbitrary T > and rewrite (39) in 11 in the integral form 

t n 

Q!l{x,t) = {Ria!^){x,t) - j "^{aijaj) {uii{T;x,t),T)dT, i<n, (41) 

ti{x,t) 
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where {Ria^){x,t) = if ti{x,t) = and 

n 

{Ria%x,t) = ^ {Pij + Kfj) {U{x,t))vf' {ti{x,t)) + Ll{ti{x,t)) (42) 

otherwise. We will use a modification of (41) where the boundary term {Ria^){x,t) does not 
depend on v°''^ . More specifically, using (41), we express vj'^ {ti{x,t)) contributing into (42) 
in an integral form. If the resulting expression for a^{x,t) still depends on v"'^, we apply 
the same procedure again. In a finite number of steps we reach the initial axis and thereby 
obtain the desired representation for aj{x,t). It is the sum of a continuous function (of Kfj, L|, 
and pij) and a finite number of integrals over characteristic curves in the EPs passing through 
cjj(r; x, t) and restricted to 11*. Let N be the maximal number of such integrals where the range 
of maximization is i < n and {x,t) G 11^. Applying the sequential approximation method, we 
easily derive the apriori estimate in 11^: 

N 

sup |af I < sup \Lf\ I sup \pij\ + sup \Kfj\ 

£,i,x,t i,t,£ „_n. \ i,j,t e,i,j,t 



i 





X exp < nT sup \aij \ I sup \pij 

y i,j,x,t \i,j,t 

This implies that a'^{x,t) is bounded uniformly in {x,t) G 11 and e > 0, as desired. 

It remains to prove that a^{x, t) ^ as e — )■ pointwise off J*. Note that, due to the support 
properties of z^ we have i;^+^'=(i) = v'^'%t) and hence Lf (t) = on \ {N^ U {t = 0}). As 
j^ei (- for all si < £2, we have the identity Lf{t) = on dU\{N^'>U{t = 0}) for an arbitrarily 
fixed £0 > and all e < Eq. Fix {x,t) € IP' \ and Eq > so that {x,t) € 11^' \ J^". Then 
a'^{x,t) is representable in the form (41) with Lf{ti{x,t)) = 0. Note that ti{x,t) € dU \ N^°; 
hence f"'^(ij(x,t)) is expressible in the form (41) again with Lj = 0. Continuing in this way, 
similarly to the above we arrive at the integral representation of where the boundary term 
depends neither on nor on v"'^. Roughly speaking, each a^{x,t) is a finite sum of integrals 
(of kind as in (41)) over characteristics in the EPs passing through LOi{T;x,t) and restricted to 
n*. We split each of the integrals into two parts, the splitting procedure being illustrated by 
example of the integral in (41): 



J dr = J dr + J dr, 

ti{x,t) [ti{x,t),t\\^eoip:,t) Aeo(a;,t) 



where A£p(a;,t) = {r G \ti{x,t),t] : {uJi{T;x,t),T) G J^°}- Since is bounded on 11 uniformly 
in £ > 0, there is C > (depending on T) such that the absolute value of the second integral is 
bounded from above by Ceq. This gives us the following apriori estimate in 11^: 

N / 



sup lafi < Ceq n'^ sup \pij \ + sup \Kf^ 

e<eo,i,x,t \i,j,t £>Oi,j,t 



X exp < nT sup \aij \ ^ n^ I sup \pi 

\ i,j,x,t \i,j,t 



sup \K, 

£>0,i,j,t 




(43) 
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Note that (43) is true with £i in place of Eq for an arbitrarily fixed £i < eo- The desired 
convergence is thereby proved. The first statement of the theorem now follows by dominated 
convergence. 

2. This part easily follows from diagonality of matrix D by the method of characteristics. 

3. Given e, denote by jB^ the continuous solution to (37) with in place of z and let H) 
denote the C (ll \ j) "-solution to (37) . Let K be an arbitrary compact and connected subset 

of n \ J. 

Claim 1. ^ w in C{K) as e — )■ 0. 

Proof of Claim. The difference w — satisfies system 

{dt + Ad^ + D) {w' -w) + F{z' - z) 
jB^{x, 0) — w{x, 0) 

wl{0,t) - Wi{0,t) 
wl{l,t)-Wi{l,t) 

Since z*^ — z in ^''(11) as e — > 0, we obtain the convergence w)'^ — > w) in £'{U.) as £ — ^ 0. It 
remains to prove that {w'^js^o converges in C (K). Note first that the solution to problem (37) 
is the sum of the solution to the version of (37) with br{x) = and g{x,t) = and the solution 
to the version of (37) with F{x, t) = 0. Since the solution to the latter problem is continuous in 
n, our task reduces to proving the convergence of {w^}e>o in C (K) in the case when br{x) = 
and g{x, t) = 0. 

Let be the union, over i < n, j < rrii, and a G [— of all EPs emanating from 
{x*j + a, 0). Fix an arbitrary i < n and eo > such that K C Il\ It is easily seen that, 
given {x,t) G K, wf(x,t) for each £ < is a finite sum of integrals over the characteristics 
in the EPs passing through uJi{T;x,t) and restricted to 11* n (see Fig. 4). Denote these by 
Cf,. . . ,Cf. Moreover, suppose that C| is given by equation x = Cj{t). 

Let fq^s denote the entries of the matrix F. Moreover, q{j) is the index of the characteristic 
curve passing through C? and s{j) is the index of the characteristic curves bounding the "tube" 
containing C|. An explicit calculation of wf{x,t) gives us equality 




where [t^ is the projection of C| onto t-axis and Pij{x,t) is a smooth function of pij. 

Furthermore, for any j < I there arc a{j) < n, /3{j) < rria, G N, and smooth 

mappings Tj^ : t — > ti^(cj(t),t), : t — > ti^{xi^,t) for 1 < p < 7, and Tj^ : t — > 
Wi^(0; Xi^,t), where Xi = ioi i < k and Xi = 1 ioi k + 1 < i < n. Geometrically speaking, for 
each j < I there is ^^^(j) ^(j) and a continuous path in containing both the curve C| and the 
point (a^*(j) /^O)'*^)' whatsoever £ < eq. Moreover, the value of 7 equals to the number of the 
"tubes" creating this path. 



= 
= 

n 
n 

i=i 
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5 



CO; (z;x,t) 




Figure 4: Construction of the set J|. 



Note that / as a function of (x, t) is constant on any connected component of 11 \ J*^". From 
now on by I we will denote the value of this function on K. Calculating of z^^j^ and changing 
variables by means of Tj^, we transform expression (45) for wf{x,t) on if to a form more 
convenient for our purpose, namely, 

where the kernels Qij{£,\x,t) are certain smooth functions. Therefore, 

I 

wl{x,t) ^(-l)^«0).«^)5^''^«^-«^'^^g,,«(^.),^(^.);x,t) as e ^ 0, 
where the convergence is uniform for all {x, t) E K. □ 

Claim 2. — it) + as e — >■ in C{K), where w and w are C(n \ J)-solutions to problems 
(37) and (38), respectively. 

Proof of Claim. On the account of Claim 1, it is sufficient to prove that — — w ^ as 
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£ — > in C{K). The difference — — w is a solution to the fohowing problem: 

{dt + Ad^ + A) {w^ -Oj^ -w) = -F{il)^ - w) 
{w^ — — w) (x, 0) = 

n 

k + 1 < i < n 

n 

iwl-wl-Wi){l,t) = ^pijit)vf-^'-^{0,t) + ri{t,v^')-n{t,v^+^), 

l<i<k. 



(47) 



Rewrite 



nit,v^ )-nit,v 
1 



iVyVi) {t, av^ + (1 - (t)^;^ +^) da ■ v"^ -^ -"^ + {nit, v^ +^) - nit, v'^+^)) 



Similarly to (39), consider an integral form of (47) on K with the boundary term not depending 
on y^^-'^^-'w^ Since — > it; as e ^ uniformly on any compact subset of IT \ J, the difference 
ri{t,v^^~^^{t)) — ri{t,v^~^'^{t)) tends to zero as £ — > uniformly on any compact subset of 
[dU \ {t = 0}) \ J. 

It remains to show that, for all i ^ j, the function 

t 

Pij{x, t)= J ^iiT'^ X, t)Jij{^, t) [w^j - iBj) (^, T)\^=oji{T;x,t) dr 

ti{x,t) 

converges to zero as e — >^ uniformly on K. Fix an arbitrary ei < £o and split up the integral 
into two parts, namely 

Pij{x,t)= J '^^+ / 

[ti{x,t),t]-Ae^{x,t) Aei(x,t) 

where A^^{x,t) = {r G [ti{x,t),t] : {uJi{T;x,t),T) G J^^}. Since Wj is uniformly bounded 
on each bounded subset of 11 and w G C(n \ J) , the absolute value of the second integral is 
bounded from above by Cei for all (x, t) E K and for some C > which docs not depend on Si. 
Furthermore, there exists £2 < £1 such that for all £ < £2 the absolute value of the first integral 
is bounded from above by £1 , due to the uniform convergence of iVj off J^^ . We hence conclude 
that for each £1 < £0 there is £2 < £1 and a constant C > such that for all {x,t) G K and 
£ < £2 we have \P[j{x,t)\ < Cei. Therefore, Qfj(x,t) — t- as £ — t- uniformly on K. □ 

To finish the proof of this part of the theorem, it remains to recall that problem (36) for w 
has a unique solution of the form w + w, where w and w are the unique solutions to problems 
(37) and (38), respectively. 

4. This part of the theorem is a direct consequence of the three preceding parts. 

□ 
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In Theorem 21 we established the existence of a delta wave solution to the problem under 
consideration. Now we identify conditions under which this solution is smoothing. In the 
following theorem conditions (i) and (ll) are supposed to be formulated correspondingly to 
problem (1), (31), (32) . 

Theorem 22 Assume that the data Aj, aij, gi, and pij are smooth functions in all their argu- 
ments and both ri and Vyri are in L°° ((0,T) x R") for any T > 0. Suppose that condition (5) 
is fulfilled. 

• Sufficiency. // condition (t) is true, then the delta wave solution to problem (1), (31), (32) 
is smoothing for any br G C[0, 1]" satisfying equalities (33). 

• Necessity. Assume that the delta wave solution to problem (1), (31), (32) is smoothing for 
any br G C[0, 1]" satisfying equalities (33). Then condition (a) is fulfilled. 

Proof. Suffi,ciency. As condition (z.) is true, there exists Tq such that .J.^, G U^" and the 

restriction of u to 11 \ 11^" equals w. The smoothingness of u will follow from Theorem 12 and 

y 

from the fact that there \sT >Tq such that w is continuous on 11 \ 11 .To prove the latter, fix 
T to be a supremum over f > such that there is x G [0, 1] and an EP jointing the point (x, t) 
with the line t = Tq. Note that u on 11 \ 11 fulfills the system of integral equations (lO)-(ll) 
where hi{t,v) = Yl'j=iPiji'l'')''^i(f ) + ^ii^^-'^^)- By the choice of T, the right hand side of (10) can 
be rewritten in the form involving integrals of kind Iijs{x,t) defined by (12). This expression 

will not depend neither on v nor on u{x, Tq). We therefore can use a transformation of Iijs{x, t) 

y 

similar to that made in (13). Since u is piecewise continuous on 11 \ 11 , the right hand side of 
(13) is a continuous function on 11 \ 11 , as desired. 

Necessity. By Theorem 21, the delta wave solution is given by the sum of a regular part w 
in C(n \ J)" and a singular part z G V'iji) which, in its turn, is the sum of strong singularities 
concentrated on the characteristic curves contributing into the set J*. By the assumption, the 
delta wave solution is smoothing for each br G C[0, 1]". This entails the smoothingness of both 
the regular and the singular parts of the solution. It remains to note that condition {u) is 
necessary for smoothingness of w. This follows by similar argument that was used to prove the 
necessity in Theorem 12. □ 
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